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On  the  Onset  of  Breakup  in  Inviscid  and  Viscous  Jets 

by 

D.  A.  Caulk  and  P.  M.  Naghdi 


Abstract.  This  paper  is  concerned  with  the  instability  of  inviscid  and  viscous 
jets  utilizing  the  basic  equations  of  the  one-dimensional  direct  theory  of  a 
fluid  jet  based  on  the  concept  of  a Cosserat  (or  a directed)  curve.  First,  a 
system  of  differential  equations  is  derived  for  small  motions  superposed  on 
uniform  flow  of  an  inviscid  straight  circular  jet  which  can  twist  along  its  axis. 
Periodic  wave  solutions  are  then  obtained  for  this  system  of  linear  equations ; 
and,  with  reference  to  a description  of  growth  in  the  unstable  mode,'  the  compari- 
son of  the  resulting  dispersion  relation  is  found  to  agree  extremely  well  with 
the  classical  (three-dimensional)  results  of  Rayleigh.  Next,  constitutive 
equations  are  obtained  for  a viscous  elliptical  jet  and  these  are  used  to 
discuss  both  the  symmetric  and  the  anti-symmetric  small  disturbances  in  the 
shape  of  the  free  surface  of  a circular  jet.  Through  a comparison  with  available 
three-dimensional  numerical  results,  the  solution  obtained  is  shown  to  be  an 
improvement  over  an  existing  approximate  solution  of  the  problem. 
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1. 


Introduction 


The  disintegration  or  breakup  of  a fluid  jet  due  to  surface  tension,  or 
so-called  capillary  instability,  has  long  been  a subject  of  interest  in  fluid 
mechanics.  We  do  not  include  here  a complete  list  of  references  on  the  subject, 
but  mention  those  that  have  some  bearing  on  the  present  paper.  The  classical 
study  for  an  inviscid  jet  is  due  to  Lord  Rayleigh  [1,2],  who  in  turn  cites  the 
static  investigations  of  Plateau  and  the  experiments  of  Savart.  Later  Rayleigh 
[3]  extended  his  earlier  work  to  include  viscosity  but,  due  to  analytical  dif- 
ficulty, explicitly  considered  only  the  case  of  negligible  inertia.  Weber  [4] 
also  examined  the  stability  of  a viscous  jet  and  sought  to  remove  the  difficulty 
encountered  in  Rayleigh's  work  [3]  by  introducing  approximations  to  the  three- 
dimensional  theory  that  take  account  of  the  'thinness'  of  the  jet.  These 
approximations  are  ad  hoc  in  nature  and  somewhat  inconsistent,  but  they  lead  to 
results  which  agree  reasonably  well  with  Rayleigh's  [1,2]  in  the  special  case 
of  an  inviscid  jet. 

In  the  references  cited  above,  the  problem  of  jet  instability  has  been 

approached  by  considering  small  perturbations  to  a uniform  cylindrical  jet  in 

* 

the  context  of  the  linearized  three-dimensional  equations.  This  procedure  leads 

to  a relatively  simple  dispersion  relation  in  the  case  of  an  inviscid  jet  [1,2], 

but  when  the  fluid  is  viscous  the  frequency  of  the  wave  motion  is  governed  by 

an  implicit  transcendental  equation  and  many  results  must  be  obtained  numerically 

(see  Chandrasekhar  [6])  or  by  approximation  [4].  Rather  than  consider  another 

such  approximation  here  we  approach  the  subject  using  a one-dimensional  theory 

of  a directed  (or  a Cosserat)  curve  in  the  form  given  by  Green,  Naghdi  and 

Wenner  [7].  The  relevance  and  applicability  of  this  approach  to  problems 

involving  fluid  jets  have  been  demonstrated  in  several  papers  by  Green  and 
* 

These  references  consider  the  temporal  instability  of  an  infinite  jet.  This 
should  be  distinguished  from  the  so-called  spatial  instability  of  a semi- 
infinite  jet  considered  by  Keller  et  al.  [5]* 


Laws  [8],  by  Green  [9,10]  and  by  Caulk  and  Naghdi  [11].  Additional  background 
on  the  theory  of  a Cosserat  curve  and  its  applications  can  be  found  in  Green 
et  al.  [7,12]. 

With  the  use  of  the  exact  linearized  three-dimensional  equations , Rayleigh 
[1,2]  derived  the  explicit  result  that  the  jet  is  unstable  only  in  the  axi- 
symmetric  mode  of  disturbance.  Inasmuch  as  the  present  development  does  not 
begin  with  the  three-dimensional  equations,  all  modes  of  disturbance  which 
occur  in  the  present  one-dimensional  direct  theory  must  be  examined  for 
stability.  It  is  because  of  this  that  in  section  2 of  the  paper  we  begin  with 
a brief  review  of  the  basic  equations  for  a straight  jet  of  elliptical  cross- 
section  which  can  twist  along  its  axis?  In  section  3,  we  derive  linearized 
field  equations  for  an  incompressible  inviscid  jet  appropriate  for  small  motions 
superposed  on  uniform  flow  of  a circular  jet.  We  show  that  the  solution  to 
these  linearized  equations  can  be  decomposed  into  two  modes,  representing  a 
symmetric  and  an  anti-symmetric  disturbance  in  the  shape  of  the  free  surface. 

The  anti-symmetric  mode  is  stable  for  all  wavelengths,  while  the  symmetric 
mode  is  found  to  be  unstable  over  a range  of  longer  wavelengths.  These  con- 
clusions are  then  compared  with  the  corresponding  exact  three-dimensional 
analysis  of  Rayleigh  [1,2]?  In  terms  of  a description  of  growth  in  the  unstable 
mode,  the  agreement  in  this  case  is  remarkably  good. 


The  constitutive  equations  for  a linear  viscous  elliptical  jet  in  the 
absence  of  twist  are  considered  in  section  4 and  are  utilized  in  section  5 


Consideration  of  a twisted  elliptical  jet  will  permit  growth  of  a general 
disturbance  which  is  not  necessarily  symmetric.  Indeed,  if  a priori 
assumptions  are  made  to  restrict  the  motion  of  the  jet,  any  conditions 
for  instability  (in  the  context  of  the  present  theory)  are  only  sufficient. 

^Here  our  results  partly  overlap  with  a recent  study  by  Bogy  [13],  who  deals 
only  with  axially  symmetric  disturbances  of  a nonrotating  jet.  In  particular, 
Bogy  considers  the  instability  of  an  incompressible  liquid  viscous  jet  of 
circular  section  and  in  the  main  his  work  is  concerned  with  spatial  instability 
of  a semi -infinite  jet  formulated  as  a boundary  value  problem. 
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2.  basic  equations. 

We  summarize  in  this  section  the  main  kinematics  and  differential  equations 
characterizing  the  motion  of  a directed  fluid  jet  in  the  form  derived  by  Caulk  and 
Naghdi  [11]  • The  jet  is  straight,  incompressible  and  homogeneous.  Recall  that 
this  characterization  of  fluid  jets  is  based  on  the  concept  of  a Cosserat  (or 
a directed)  curve,  hereafter  designated  as  ft.  Such  a one-dimensional  directed 
medium  comprises  a material  line  and  a pair  of  directors  attached  to  every 
point  of  the  material  line. 

Let  the  particles  of  the  material  line  of  ft  be  identified  with  the  convected 

coordinate  let  c,  the  curve  occupied  by  the  material  line  of  ft  in  the  present 

configuration  at  time  t,  be  described  by  its  position  vector  r relative  to  a 

fixed  origin;  and  let  d (or  = l,2)  stand  for  the  pair  of  directors  at  r.  Then, 

a ~ 

a motion  of  the  directed  curve  R is  specified  by 


r.r(5,t)  , 4,-d^ft.t)  . 

The  velocity  and  the  director  velocities  are  defined  by 


(1) 


v = r , w = d 

~ ~ ~oi 


(2) 


where  a superposed  dot  designates  the  material  time  derivative  holding  5 fixed. 

For  the  purpose  of  displaying  the  details  of  the  kinematics  of  a straight 
jet,  including  the  rotation  of  the  directors  in  a plane  normal  to  the  jet  axis, 
we  introduce  a fixed 'system  of  rectangular  Cartesian  coordinates  (x,y,z)  with 
the  z-axis  parallel  tc  the  jet.  Further,  let  the  unit  base  vectors  of  the 
rectangular  Cartesian  axes  be  denoted  by  (i,j,k)  and  introduce,  for  later 

rsw  /-w  rw 

convenience,  the  additional  base  vectors 

e = i cos  9 + j sin  0 , e0=  - i sin  0 + j cos  9 , e_  = k , (3) 

where  9 is  a smooth  function  of  z and  t.  We  assume  that  the  directors  are  so 


restricted  that  they  describe  an  elliptical  cross-section  of  smoothly  varying 
orientation  along  the  length  of  the  jet  and  that  at  each  z = const.,  the  base 
vectors  e^  and  e^  lie  along  the  major  and  minor  axes  of  the  ellipse,  respectively. 
Then,  the  angle  0,  called  the  sectional  orientation,  specifies  the  orientation 
of  the  cross-section  as  a function  of  position.  With  this  background,  henceforth 
we  restrict  motions  of  the  directed  curve  R such  that  in  the  present  configuration 
at  time  t, 

r»z(5>tk3  , ^-0^  , (U) 

where  01  and  measure  the  semiaxes  of  the  elliptical  cross-section.  The 
velocity,  the  acceleration,  the  director  velocities  and  the  director  accelera- 
tions assume  the  form^ 


v = ve  , v = z , v = ve.  , 

^3  rsj 

Si  = 0i(ci£i  + u’iS2)  ' Ss  = 


where 


51  = (^1  + ^1  " u,l)0l2l  + (2Ci»i  + «i  )0i%>  > 

52  = ( ^1  + 4 ' 4 )02^2  " ^^2  + *2  )02~l  ’ 


01  01^1  * 02  02^£ 


and  and  ^ represent  the  rotational  components  of  the  director  velocities  in 
the  plane  normal  to  the  jet  axis. 

The  condition  expressing  the  incompressibility  of  the  fluid  medium  is 
given  by 


W«2  * ° 


It  should  be  noted  that  (6)  and  (7)  represent  values  for  w and  w in  the 
present  configuration  and  are  not  obtained  by  direct  differentiation^?!  the  cor- 
responding present  values  (4)p  _ for  d ; see,  in  this  connection.  Caulk  and 

Naghdi  [11]. 
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and  the  differential  equations  of  motion  for  the  elliptical  jet  described  above 


are  (see  Caulk  and  Naghdi  [11]): 


A 

dn 


i?  = tpz  +020lZh(0l’8>2>9z)  +^2zh(02J0i.9z)  + 9z(02-01)m(?>1,02,0z) 


+ np*0102v}eg  , 


(10) 


dp 


A1 


(-§—  ■ = c-p-0i02h^i^2’0z^  + ^TTp^i02^i  + ci"oui^ 


Si 


+ [0^m(01,02,0z)  + £ + 2£1u)1)]e2  , (11) 


*2 

dp 


" 22/z')02  = t-P  '^102h(02,0l’ez^  + '**~  npX0201^2  + £2  "“>2^% 


2 ^ o # 

- {02m(02,01,ez)  +i  np  0^01(U)2 + 2C2-jj2)} 


Si  > 


dd 


, A V Afy 

z e_  x n + d xu  + -rz-  X P =0 

~3  ~ ~0f  ~ dc,  ~ 


(12) 

(13) 


where  z'  = 3z/d§  and  a subscripted  z denotes  d/dz.  The  kinetic  quantities  n, 


tt  and  p are  specified  by  constitutive  equations , p is  related  to  an  average 

/v»  rsj 


* . 


pressure  over  the  cross-section  of  the  jet  and  p is  the  three-dimensional 
density  of  the  fluid.  In  addition,  h and  m arise  due  to  the  constant  surface 
tension  T and  are  given  by 


h(0 


1 2 z' 


m 


z.' 


r2^  2 

1 q cos  x IX  , 

(14) 

°0 

.2tt 

j q sin  2x  dx  , 

(15) 

where 
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j 


q = 


-I 

i (I  G^sin  X cos  X ) ~ (^l^ozs^n * 020^zcor.2x)  J2  1 J^sin^X  + 02cos2x} 

X { (0^sin2x  + ^c°s2x)  (<^1zz02cos2x  + ’ [ 6z(02_;Sl j !zSin  X COS  X 

" 0l029z^  “ 2^0102z  "020lz^Sinx  COSX  " 9Z(02COSX  + ^isin2x)J 

X [(0202z-010lz)sinXcosx+  01020z] 

- 0102[(0lzcosx-  02ezsinx  )2+  (02zsinx+  010zcosx  )2  + l]}  , (16) 

and  for  convenience  we  have  let  * stand  for  the  triple  U ,#>  ,&  ) in  the 

a a az  azz 

arguments  of  h and  m.  Also,  we  note  the  fact  that  m satisfies  the  condition 

m(01,02,O)  = 0 . (17) 

To  complete  the  above  system  of  equations,  we  must  add  an  expression  for 
6>  or  the  rate  of  sectional  rotation  of  the  jet.  When  the  jet  is  non- 


circular,  this  relation  is 


6 = 


2 2 
V2  ~ Vl 
*2  *2 
02'01 


(18) 


For  later  reference,  we  record  here  the  expression  for  the  mechanical  power 
F per  unit  mass  of  an  incompressible  jet,  namely 


A Sv  a a 3w 

1 P = n • T7  + tt3  • v + pa  • , 


(19) 


wnere 


\ = z'p  . 


(20) 


At  a point  where  the  jet  is  circular  (01  = 02),  there  is  no  preferred  orientation 
of  the  cross-section. 
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3.  Small  motions  superposed  on  a uniform  flow  of  nn  inviscid  .jet 

First  we  consider  an  inviscid  jet  characterized  by  the  constitutive 
assumptions 


A /Vv  A-v 

n = 0 , tt  = 0 , p = 0 


In  this  case,  (10)  to  (12)  reduce  to  the  five  scalar  equations 


(21) 


"P  ^i^2^Vt  + VVJ  = - pz  “^20lzh^l,02,0z^  ‘ ^l42zh^2,0l,0z^  ’ ez^2'Vm^l,?2’9z^  * 


-3  2 2 

i "P  ^2^lt  + vClz  + tl-u,l^  = P + ^n0o’n(01  ,0O,6 


(22) 


lv2“v’'l,w2’uz/  ’ 

1 ^ 

4 up  0102^lJult  + Vtulz  + ^®1U)1^  =~m^x’32’®"^  ’ 

i TTP^25>1(t2t  + V^2z  + ^ ■'i,2)  = P + di02h(02»iir0z)  » 

u up  ^1^2 ^2t  + V(U2z  + 2=2(u2'  =“m^2,^l’®z^  ’ 

where  we  have  expressed  all  functions  in  terms  of  the  current  position  z of 
the  material  particle  | at  time  t.  The  set  (22)  is  completed  by  adding  the 
incompressibility  condition  (9)  together  with  the  kinematic  relations  (8)^  ^ 
and  (18).  Consider  now  the  simple  solution 

v = vq  , 01  = 02  = a , Cx  = »2  = 0 > u>i  = uu2  = 0 , p = u&T  , (23) 

which  satisfies  exactly  the  above  system  of  differential  equations  and  where 
vq  and  a are  constants.  This  solution  represents  a uniform  circular  jet 
moving  with  constant  velocity  v . With  a suitable  choice  of  reference  frame 
we  may  let  v =0;  then,  due  to  the  rotational  symmetry  of  the  directors  in 


the  plane  of  the  cross-section,  we  may  (at  least  for  this  solution)  take 


*z  ■ 0 ■ 


t3**> 


without  loss  in  generality. 


8. 


We  now  examine  small  motions  superposed  on  the  uniform  flow  represented  by 
the  solution  (23) • This  naturally  leads  to  a discussion  of  jet  instability  and 
breakup  (or  disintegration)  as  each  is  generally  understood  in  the  literature. 
Accordingly,  we  shall  determine  a linearized  version  of  the  governing  equations 
(22),  (8),  (9)  and  (l8)  appropriate  for  small  deviations  from  the  motion  (23). 
We  proceed  in  this  manner  owing  to  its  relative  simplicity  and  wide  use,  but 
note  that  any  results  arising  from  such  a treatment  shall  be  necessarily 
restricted  by  the  limited  scope  of  a linearized  stability  analysis. 

Consider  small  deviations  from  the  motion  (23)  in  the  form 


0 


a 


a + ?* 
a 


v 


v + v 
o 


p = rraT  + p 


(25) 


and  retain  only  linear  terms  in  quantities  represented  by  symbols  with  superposed 

tildas  ("~")  in  all  equations.  Drawing  upon  the  discussion  preceding  (24) 

we  take  v =0  and  0 = 0 in  the  unperturbed  flow.  Then,  in  keeping  with  the 
o z 

linearized  procedure,  we  assume  that  0 is  small  in  the  perturbed  flow  and 

z 

neglect  its  square  and  product  with  quantities  having  tildas.  If  we  linearize 
(16)  in  this  manner,  with  the  help  of  (14)  and  (15)  we  obtain 


ttT  ttT  r ^ 2.**  1 2~  1 ^ ^ _ 

h ="T  + ~2  C4  a 0lzz  + 4 a 02zz  + £ ^2"0l)+02]  ’ m = 0 • 


(26) 


Introducing  (25)  and  (26)  into  (22),  (8)  and  (9),  neglecting  squares  and 
products  of  small  quantities  and  then  dropping  the  tildas  for  simplicity,  we 
are  left  with^ 


^Since  neither  0 nor  its  derivatives  appear  in  any  of  the  other  linearized 
equations  (27) > we  do  not  record  a linearized  counterpart  to  (18). 


rrp  a vt  = - Pz +rrT(01  + 0g)z  , 


i np  a^Clt  = P + TTT[^-  (30-j^  + 02)zz  - 02  +'i(0?  - 0-l)]  > 


4 nP  a^^2t  = p +ttT^V  ^2  + ^l\z  “ ^2  + ^1  “ $2^  ’ 


£l  + £2  + Vz  = ° 3 


a^l  ^lt  ’ a^2  02t  5 


“it  = ° ’ ^t  = ° • 


In  order  to  examine  solutions  to  the  system  of  partial  differential  equa- 
tions (27) » it  is  convenient  to  introduce  the  change  of  variables 


0=i(01  + 02)  > 6 = |(01-02) 


Substituting  (28)  into  (27)  and  then  adding  (27)2  and  (27 5 with  the  help  of 


(27)c 


>6,7,8 


we  obtain 


* 2 

p tte  v,  = - p + 2ttT0  , 
t z z 


^p*rra30tt  = Pt  +nT[a20zz  - 0]  , 


20,  + av  =0 
T>  2 

Alternatively,  if  we  use  (28)  in  (27)  and  subtract  (27  )^  from  (27 )g,  we  have 

6tt  ■ -§%  <«28M  - 36)  • 

P a 

Hence,  in  terms  of  the  variables  (28),  the  set  (27)  decouples  into  (29)  and 
(30)  and  we  can  find  solutions  to  each  separately.  First,  we  consider  (29). 
Elimination  of  p and  v among  these  equations  yields 

^ttzz  " 20tt  = ~ (*zz  + a ^zzzz’ 
p a 


• |*^%**S^ 


as  a differential  equation  in  0 only.  We  examine  solutions  of  (3l)  in  the 
form 


0(z,t)  = fQ  (ko)exp[i(aQt  - kQz)]  , 


(32) 


from  which  follows  the  dispersion  relation 

2 4T 


ao  = 


.2  2 . 

( °a  )i-2a2 

pV  Wa2'  0 


(33) 


From  (33)  it  is  clear  that  the  wave  motion  (32)  is  unstable  for  wave  numbers 
satisfying 


•V>2  < 1 • (34) 

Returning  to  (30),  we  consider  periodic  solutions  to  this  equation  in  the 
form 


6(z,t)  = f2(k2)exp[i(c2t  - k2z)]  . 


(35) 


The  resulting  dispersion  relation  is 


= (k2a2  + 3) 

P a 


2 

32 


Now,  from  (28)  we  have  for  0 and  02 


01  = 0 + 6 , 02  = 4 - 6 , 


(36) 


(37) 


so  that  the  general  motions  consists  of  two  parts:  When  6=0,  (37)  gives 
01 = 02  = 0 and  the  circular  cross-section  of  the  jet  remains  circular  in  the 
perturbation.  We  call  0 the  symmetric  mode.  When  0 = 0,  we  have  0^  = -02  = 6 
and  the  cross -section  of  the  perturbed  jet  is  an  oscillating  ellipse, 
exchanging  alternatively  its  major  and  minor  axes.  This  we  call  the  anti- 
symmetric mode. 


In  summary,  a small  disturbance  to  the  motion  (23)  can  be  decomposed 


! 

I 


into  two  modes  of  vibration:  (i)  a symmetric  mode  whose  frequency  is  governed 

by  (33)  and  (ii)  an  anti -symmetric  mode  whose  frequency  satisfies  (36).  The 

latter  mode  is  stable  for  all  wavelengths,  while  the  symmetric  mode  is  unstable 

for  wavelengths  satisfying  (3*0-  It  is  worth  noting  here  that  if  the  rotations 

(u  and  the  twist  9 are  set  equal  to  zero  at  the  outset  in  the  equations  (10) 
a z 

to  (13 ),  then  the  results  and  conclusions  of  this  section  would  remain 
unchanged. 

Before  closing  this  section,  we  make  a comparison  with  the  results  of 
Rayleigh  [1,2]  who  examined  small  deviations  from  the  uniform  flow  oi  an 


inviscid,  incompressible  fluid  in  a straight  circular  jet.  For  this  purpose, 
we  introduce  cylindrical  polar  coordinates  (r,X,z)  such  that  the  z-axis  lies 
along  the  axis  of  the  jet.  Rayleigh  considered  disturbances  from  a circular 
jet  of  radius  a in  which  the  free  surface  had  the  modal  forms 


r = a + b(z ,t )cos  n X , 


where  n is  an  integer  and  b is  small.  With  (38)  as  his  basic  assumption, 
Rayleigh  examined  solutions  to  the  linearized  three-dimensional  equations  for 
which  b(z,t)  has  the  form  (32).  The  resulting  dispersion  relations  for  each 


integer  n (as  recorded  in  Lamb  [lU,  §274])  are 

_ _ l'(k  a ) [k2a2  + n2  - l]k  a 

2 T nv  n n n 


a = r 

n *3 


I (k  a, 
n n 


where  I is  the  modified  Bessel  function  of  order  n.  The  only  value  of  n that 
n 

can  lead  to  unstable  wave  motion  is  n = 0,  corresponding  to  an  axially  symmetric 
disturbance.  With  n = 0,  (39)  becomes 

_ I ' (k  a)[k2a2  - l]k  a 

2 T 000  o /\ 


°o  *3 
P 3 


I (k  a ’ 
0 0 


indicating  a range  of  unstable  wavelengths  consistent  with  (34).  The  unbounded 
growth  of  disturbances  in  this  range  leads  to  an  eventual  disintegration  of  the 
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jet.  For  purposes  of  comparison,  we  tabulate  in  Table  1 values  for  a from 
(33)  and  (40)  over  the  entire  range  of  unstable  wavelengths.  Of  the  wave- 
lengths in  this  range,  the  one  that  corresponds  to  the  greatest  magnitude  of 
aQ,  and  hence  the  most  rapid  rate  of  growth  in  the  disturbance,  will  tend  to 
dominate  the  disintegration  process.  The  value  obtained  by  Rayleigh  for  this 
wavelength  corresponds  to 

k2a2  = 0.4858  , (4l) 

whereas  from  (33)  we  find 

k2a2  = 0.4853  (42) 

o 


yields  maximum  growth. 

Setting  n = 2 in  (39),  we  obtain  the  counterpart  to  (36)  in  the  linearized 
three-dimensional  theory,  namely 


2 T l2^k2a^k2&2  + 3'lk2a 
a2  = 

P 8 


I2(k2a) 


(43) 


In  line  with  (36),  the  motion  governed  by  (43)  is  stable  for  all  wavelengths. 

$ 


2 2 

Expanding  (43)  in  powers  of  k a yields 


°2  = ”*^3  (k2&2  + 3 ) [1  + Xo  k2&2  + ' * * 1 • 

P a 


(44) 


A comparison  of  this  relation  with  (36)  suggests  the  appropriateness  of  the 
latter  for  fairly  long  wavelengths. 
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(— 2_)-i 

v ' W t 


koa 

Directed  Jet 

j Rayleigh 

i 

0.00 

1 

0.0000 

0.0000 

0.03 

0.1536 

j 0.1536 

1 

0.10 

0.2107 

0.2108 

0.20 

0.2793 

0.2794 

0.30 

0.3181 

0.3182 

0.40 

0.3381 

0.3332 

0.50 

0.3429 

0.3432 

0.60 

0.3341 

0.3344 

0.70 

0.3107 

0.3111 

0.80 

0.2696 

0.2701 

0.90 

0.2010 

0.2015 

1.00 

0.0000 

0 . 0000 

Table  1:  Coir.parisori  of  the  frequency  a0  (in  non-dimensional 
form)  over  the  range  (3*0  of  unstable  wavelengths 
as  predicted  by  Eq.  (33)  of  the  direct  theory  of 
jets  and  by  Eq.  (40)  due  to  Rayleigh  [1,2]. 


-*  * * »-~Y*  V' 


A viscous  elliptical  .iet  without  rotation  or  twist. 


In  this  section  we  consider  a jet  of  an  incompressible  linear  viscous 

fluid;  and,  in  view  of  the  results  for  the  stability  of  an  inviscid  jet,  we 

limit  the  discussion  to  motions  in  which 

u>l  = m2  = 0 > 0 = const.  ( 

To  account  for  the  viscosity  of  the  fluid  medium,  we  must  provide  appropriate 

A ^\v 

constitutive  equations  for  the  quantities  n,  n and  p . Much  of  the  develop- 

n>J  r*0  r* 

ment  of  this  section  is  similar  to  that  of  Green  [10]  in  which  the  jet  is 
restricted  to  be  circular. 

Referred  to  the  orthonormal  basis  (e, ,eoJe_)  introduced  in  section  2,  the 

'VL 

response  functions  can  be  expressed  in  terms  of  their  components  in  the  form 

A Ai  Aq , Aori  Aori  / 

n = n e.  , n = tt  e,  , p = p e.  • ' 

/yl  r+J 

As  a consequence  of  the  symmetry  of  the  assumed  flow,  (10)  to  (12)  suggest 


that  we  put 


Ai  Ae 

n = n = 0 , 


A12  a13  *21'  ^3 
TT  =TT  =TT  =TT  =0  , 


A12  ai3  Aei  ^3  „ 

p = p = p = p = o 


With  the  help  of  (4)  to  (6)  and  (46),  the  mechanical  power  (19)  reduces  to 


X D a3  **11,.  , ^*22.  .£*11,.  A' 
— P = n vz  + tt  C1  + tt  C2  + P £lz  + p 


®2z  ’ 


where  we  have  set 


A*oar  . ,\xat/  A*aar  _ . /Va  __  \ 

tt  =0VTT/zj+0p  jP  — 0 P \t\o  sum  on  ot)  • 
a Qfz  or 


We  seek  to  characterize  the  linear  viscous  property  of  the  fluid  in  appropriate 
constitutive  equations  for  the  one-dimensional  functions 


v ■ ’T*  « t v *>%  r — 

ii.i i a 


(50) 


A3  Axil  Ax22  Axil  ^*22 
n ,n  ,n  ,p  ,p 


To  this  end,  we  assume  that  the  quantities  (50)  are  linear  functions  of  degree 
one  in  the  kinematic  variables 


V^L’^’^lz’^z 


(51) 


with  coefficients  that  depend  upon  and  0^.  Hence,  we  take 


nA--xfS+xfSz  , 


1 (no  sum  on  a) 


where  are  functions  of  01  and  0^  and  we  have  used  (9)  to 

A*(V(y  A#/y/y 

eliminate  v from  the  expressions  for  n and  p 
z 

A three-dimensional  linear  viscous  fluid  is  isotropic.  In  order  that  the  one- 
dimensional  theory  under  consideration  reflect  the  symmetry  properties  of  the  fluid 
and  the  geometry  of  the  jet,  we  impose  the  requirement  that  under  the  transformations 


(52) 


z - - z 


V -*  - V 


, C £ , <f>  - t , 

a a a 


(53) 


the  mechanical  power  (48)  remain  invariant.  Consequently,  under  (53)  the 
functions  (50)  must  transform  according  to 

A3  . A3  £*aa  A *aa  £*aa  Ax-ca 
n-*n  , rr  -*tt  jP  - P 


(54) 


Hence  the  relations  (52)  reduce  to 
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■ — 


r 


* A • 4 v ^ * 


(55) 


I , 

r 


A- 


n3  = \ v\  + , 

I X 2^of 


TT*a“  = (n°  SUm  °n  “)  > 


A* 


Ora 


= x“^Cfl  (no  sum  on  or)  . 


7 b0z 

In  order  to  determine  explicit  values  for  the  coefficients  in  (55)j  we 

recall  briefly  some  aspects  of  an  approximation  procedure  for  rod-like  bodies 

in  the  three-dimensional  theory.  A detailed  development  of  this  procedure  can 

be  found  in  Green  et  al.  [12]  and  a brief  outline  is  included  as  an  appendix  to 

Caulk  end  Naghdi  [ll].  Without  going  into  detail,  we  recall  that  the  developments 

in  Green  et  al.  [12]  are  based  on  an  approximation  for  the  position  vector  to  the 

* ... 

material  points  in  the  rod-like  body  and  involve  integration  of  the  three-dimensional 

equations  through  its  cross-section.  Let  the  material  points  be  identified  with  the 

i 3 

convected  coordinates  9 (i=  1,2,3)  and,  for  convenience,  set  9 = £.  Further,  let 

p denote  the  position  vector  of  a typical  point  at  time  t.  Then, 


op 


p = p(ea,§,t)  , g , = —T  , 

~ ~ ~1  U01 


(56) 


.i  < s S.  • g.  > g * g*  - 
ij  ~i  ~ ~0  0 


i in  : 

, g = g 


, g = det  gi>.  , 


where  g.  and  g^  are  the  covariant  and  contravariant  base  vectors,  g. . is  the 


metric  tensor,  its  inverse  and  6^  is  the  rCronecker  delta.  In  the  present 

J 


context  the  fluid  is  assumed  to  occupy  a region  of  space  in  the  neighborhood 
of  the  curve  9Q=  0,  bounded  by  the  free  (material)  surface 


/ nl\2  , ' , 

( 9 J + ( 9 ) =1  5 


(57) 


where  we  identify  0a  = O with  the  z-axis.  With  the  help  of  (4), 


It  should  be  mentioned  that  the  one-dimensional  equations  that  result  from 
this  procedure  can  be  brought  into  1-1  correspondence  with  the  theory  of  a 
directed  curve. 
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the  approximation  for  the  position  vector  mentioned  above  leads  to 


P = z(a>t)e3  + + e202(S,t)r2  . (58 

It  follows  from  an  examination  of  the  second  and  third  terms  of  (58)  that 
(57)  represents  an  elliptical  cross-section  with  semiaxes  3 . The  velocity  v*, 
which  is  the  material  time  derivative  of  p,  is  given  by 

{ = v£3  + 6^i^i + ’ (59: 

where  we  have  used  (4),  (8)  and  (45).  For  an  incompressible  linear 
viscous  fluid,  the  determinate  part  of  the  stress  response  is  given  by 

= g2TlJg,  , t = rfv*.  • g.  + v*  • g ] , (60; 

where  t.  . and  t1J  are  the  covariant  and  contravariant  components  of  the  stress 
J 

tensor,  n is  the  shear  viscosity  and  a comma  denotes  partial  differentiation 
with  respect  to  Q1.  Jve  may  now  use  (57),  (59)  and  (60)  in  the  usual 
expressions'  for  the  quantities  (46)  in  terms  of  integrals  over  a cross- 
section  of  the  rod-like  body.  The  results  of  this  rather  long  but  routine 
calculation  are 


n = 2y.rr  01i2vz  , 


tt  J = 2un  a.$or  , 


~ * un  *1*2*q&xz  ’ 


| (no  sum  on  a) 


so  that  we  may  identify 


See,  for  example,  equation  (A13)  of  Caulk  and  Naghdi  [11]. 


\ ^ = 2^rj  0^&2  > = 0 , ~ 2^,77  0^2  9 

11  _ x 3 ,22  _ i ^ 3 ,12  __  ,21  _ ,12  _ ,21  _ 

Krj  — 4 ^2^1  * A.y  — 4 M'TT  0 ^02  * ^4  ~”  ^4  — ^y  ” ^y  ” ^ 

in  (55). 

Adopting  the  values  (62)  and  using  the  constitutive  equations  (55)  in 
tiie  field  equations  (10)  to  (13),  we  obtain  the  governing  differential  equa 
tions  for  a linear  viscous  jet,  namely 


* 

tip  0102(vt  + vvz)  = - Pz  - 020lzh(01,02)  - i102zh(02.01) +2p,n(?>102vz)2 


4 np  0i02(Clt  + vClz  + Cj)  +2htt  0^2^  = P + 0102h(01,02)  + £ 


4 np*0^01(C2t  + vC2z  + ;|)  + 2^rr  0X02C2  = P + 0102h(02 ,*,)  + £ nn(^01C2z)a 


where  use  has  been  made  of  (45)  and  (17)  and  we  have  let 


h(0^,02)  - h(0^»02,O) 

The  set  (63)  is  completed  by  adding  (8)  and  (9)*  Apart  from  differences 

in  notation,  we  note  that  (63)-,  o o reduce  to  those  given  by  Green  [10, 

,0 

Ics.  (6.3)  and  (6.4)]  in  the  special  case  of  a circular  cross-section  (01 = 02) 


and  in  the  absence  of  gravity. 


5.  Small  motions  superposed  on  a uniform  flow  of  a viscous  jet. 

The  motion  (23)  satisfies  the  differential  equations  (63),  (8)  and  (9); 
and,  hence,  it  also  represents  an  exact  solution  for  the  viscous  jet  discussed 
in  section  4.  Since  in  the  case  of  an  inviscid  jet  (see  section  3)  the 
superposed  rotations  uT^  had  no  effect  on  the  resulting  differentiation  equa- 
tions, for  simplicity  we  assume  that  uT  = 0 here  and  consider  small  motions 

a 

superposed  on  the  uniform  flow  (23)  in  the  form  specified  by  the  first  two 
of  (25)  and  the  fourth  of  (25)-  In  a manner  similar  to  that  employed  in 
section  3>  we  again  neglect  squares  and  products  of  quantities  represented  by 
symbols  with  superposed  tildas  in  (63),  (8)  and  (9)-  After  setting  vq = 0 
without  loss  in  generality,  the  resulting  linearized  equations  are 


^P*a2vt  = - Pz  +nT(01  +02)z  + 2^a2vzz  , 

2 

t-  TTp"a4$lt  + 2-^a2^  = p + ttT[^-  (3?1  + )M  - S,  + £(%  Clzz  » 


i no+aUC2t  + 2^a2;2  = P+nT[^-  (3*2  + ^)2Z  - *2  + ^ - ig ) 3 + £ 5 


2. 


(65) 


v„  + Ci  + ;2  - ° , 0lt  aCx  > 02t  a*2  * 

where  again  for  convenience  the  tildas  have  been  omitted. 

Again  we  utilize  the  change  of  variables  (28)  and  by  adding  (65)2  and 
(65 ) ^ ’ with  the  help  of  (65)^  ^ we  obtain 

*2  2 
' np  a vt  = - pz  + 2rrT^z  + 2p,na  , 

V np  a°itt  + 2p,naZ>t  = p +rrT[a2^5,5,  - jJ  + i > 


zz 


23,  + av,  = 0 
"t>  z 


Subtraction  of  (65)^  frojn  (65)2*  after  using  (28)  and  (65)^£>  yiexd< 
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(66) 


(67) 


2 P a36tt  + ^a6t  = T(a26zz  ■ 36)  + \ p.a^6tzz  • 

Thus,  the  linearized  system  (65)  is  decoupled  through  the  change  of  variables 
(28),  just  as  in  the  case  of  the  inviscid  jet.  Again  the  solution  of  (65) 
will  have  the  form  (37)  and  decompose  into  a symmetric  mode  0 and  an 
anti-symmetric  mode  6-  Eliminating  p and  v among  (66)^  g y ^or  the  symmetric 
mode  we  obtain 


u a20,,  - 20,.  = -%  (4  a20  - 60  ) +-7-  (a20  +0  ) 

4 ttzz  tt  * 4 zzzz  zz  t * zzzz  zz 

P Pa 


(68) 


For  solutions  of  (68)  to  be  of  the  form  (32 ),  aQ  and  kQ  must  satisfy 


2 * '8  o 
/ . \2  a 0 

(iao}  1,22 


/l  , 2 2 _ \ 2 2 

(«  kna  +3)koa 


8 V +1 


* 8 

(.  \ _ 2p  a 

(iao}  " 1 ,22 


T ,,  , 2 2x.  2 2 

(1  - k a )k  a 
' 0 0 


5 "oa  +1 


(69) 


It  follows  from  (69)  that  iOQ  can  be  real  and  positive  if  and  only  if 


,2  2 , , 
k a <1 
o 


(70) 


Hence  the  range  of  unstable  wavelengths  for  the  symmetric  mode  is  precisely  the 
same  as  for  an  inviscid  jet.  The  effect  of  viscosity  in  the  present  case  is  to 
diminish  the  magnitude  of  aQ  over  the  range  (70)  and  therefore  retard  the 
impending  disintegration  of  the  jet.  As  in  section  3>we  consider  solutions  for 
the  anti-symmetric  mode  in  the  form  (35)-  Using  (35)  in  (67)  we  obtain 
the  relation 

(io  )2+-£g  (k2a2+8)^ia2) +~^3  <4^3)  = 0 (71 ) 

pa  p aJ 

between  and  kg.  From  (71)  one  can  show  that  i Og  has  8 negative  real 
part  for  all  values  of  kg.  This  indicates,  as  in  the  case  of  the  inviscid  jet, 
that  the  anti -symmetric  mode  is  stable  for  disturbances  of  all  wavelengths. 

The  effect  of  viscosity,  however,  is  to  damp  the  disturbance  in  proportion 
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to  the  magnitude  of  the  shear  viscosity  p,.  Critical  damping  for  a given  value 
of  kg  corresponds  to 


A 


2 


Z?  = * 3) 

P*Ta  (k2a2  + 8)2 


(72) 

* 


where  we  have  introduced  the  non-dimensional  parameter  A for  later  convenience. 

We  close  with  a comparison  of  certain  results  of  this  section  with  those  of 
a similar  investigation  by  Weber  [^]  who  has  examined  small  axially  symmetric 
perturbations  to  uniform  flow  of  a cylindrical  viscous  jet  using  an  approximate 
form  of  the  linearized  Navier-Stokes  equations.  Weber's  procedure  employs 
specific  assumptions  on  the  variation  of  stress  and  velocity  in  the  cross- 
section  of  the  jet  and  ignores  all  but  the  axial  component  of  momentum.  This 
approach  leads  to  a one-dimensional  reduction  of  the  three-dimensional  equations 
and  corresponding  to  (69)  gives 


(iao)2  ♦ S$2(i»0)  - -fj 

pa  2p  a 


,,  ,2  2 2 2 
(1  - k a )k  a 
v 0 0 


(73) 


A plot  of  iaQ  versus  kQa  is  given  in  Fig.  1 for  various  values  of  the  parameter 
A over  the  range  (70)  of  unstable  wavelengths,  using  both  (69)  and  (73)-  It 

can  be  seen  from  this  graph  that  the  difference  between  the  results  of  each 
approach  is  greatest  for  an  inviscid  jet  (A  = 0)  and  gradually  diminishes  with 
increasing  viscosity,  other  things  being  equal. 

On  the  basis  of  the  close  agreement  (Table  l)  with  the  exact  three- 
dimensional  analysis  of  Rayleigh  [1,2]  for  an  inviscid  jet  (A  = 0),  it  is 
reasonable  to  infer  from  Fig.  1 that  for  a viscous  jet  the  results  of  this 
section  constitute  an  improvement  over  the  approximate  treatment  of  Weber  [4]. 

In  support  of  this  inference,  we  appeal  to  some  numerical  results  recorded  in 


¥■  2 

The  parameter  A can  be  recognized  as  twice  the  ratio  of  the  Weber  number  to 
the  square  of  the  Reynold's  number. 
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Chandrasekhar  [6]  which  are  based  on  the  implicit  dispersion  relation  obtained 
from  the  linearized  Navier-Stokes  equations t Figure  2 shows  a magnified  portion 
of  one  of  the  curve  pairs  in  Fig.  1 corresponding  to  A =0.5  along  with  points 
obtained  from  the  tables  in  Chandrasekhar  [6]  . We  show  the  region  near 
maximum  aQ  in  view  of  its  importance  in  the  breakup  process.  The  fact  that  the 
theory  of  a directed  fluid  jet  offers  an  improvement  over  Weber's  results  is 
clearly  evident. 


Acknowledgment . The  results  reported  here  were  obtained  in  the  course  of  research 
supported  by  the  U.S.  Office  of  Naval  Research  under  Contract  N00014-76-C-0474, 
Project  NR  062-531+,  with  the  University  of  California,  Berkeley. 


References 


Lord  Rayleigh,  On  the  Instability  of  Jets.  Proc.  Lond.  Math.  Soc.  10, 

4 (1879). 


2. 

3. 

4. 

5. 


8. 

9- 

10. 

11. 

12. 

13. 

14. 


Lord  Rayleigh,  On  the  Capillary  Phenomena  of  Jets.  Proc.  R.  foe.  29, 

71  (1879). 

Lord  Rayleigh,  On  the  Instability  of  a Cylinder  of  Viscous  Liquid  Under 
Capillary  Force.  Phil.  Mag.  Ser.  5 145  (1892). 

C.  Weber,  Zum  Zerfall  eines  Flussigkeitsstrahles . ZAMM  11,  136  (1931). 

J.  B.  Keller,  S.  I.  Rubinow  and  Y.  0.  Tu,  Spatial  Instability  of  a Jet. 
Phys.  Fluids  16,  2052  (1973). 

S.  Chandrasekhar,  Hydrodynamic  and  Hydroma gnetic  Stability.  Oxford: 
Clarendon  Press  (1961). 

A.  E.  Green,  P.  M.  Naghdi  and  M.  L.  Wenner,  On  the  Theory  of  Rods. 

II.  Developments  by  Direct  Approach.  Proc.  R.  Soc.  Lond.  A337,  485 
(1974). 

A.  E.  Green  and  N.  Laws,  Ideal  Fluid  Jets.  Int.  J.  Engng.  Sci.  6,  317 
(1968). 


A.  E.  Green,  Compressible  Fluid  Jets. 
189  (1975). 


Arch.  Rational  Mech.  Anal.  59, 

Int.  J.  Engng.  Sci. 


A.  E.  Green,  On  the  Nonlinear  Behavior  of  Fluid  Jets 
14,  49  (1976). 

D.  A.  Caulk  and  P.  M.  Naghdi,  The  Influence  of  Twist  on  the  Motion  of 

Straight  Elliptical  Jet.  To  appear  in  Arch.  Rational  Mech.  Anal.  = 
Report  No.  UCB/AM-77-5>  Office  of  Naval  Research  Contract  N00014- 
76-C-0474,  Univ.  of  Calif.,  Berkeley  (June  1977). 

A.  E.  Green,  P.  M.  Naghdi  and  M.  L.  Wenner,  On  the  Theory  of  Rods. 

I.  Derivations  from  the  Three-Dimensional  Equations.  Proc.  R.  Soc. 
Lond.  A337,  451  (1974). 

D.  B.  Bogy,  Use  of  One-Dimensional  Cosserat  Theory  to  Study  Instability 
in  a Viscous  Liquid  Jet.  Phys . Fluids  21,  190  (1978). 

H.  Lamb,  Hydrodynamics , 6th  edn.  Cambridge  University  Press  (1932). 


24. 


T * • 4 "*•'7  A.  ' 


UNCLASP IFIED 


SECURIT>  CLASSIFICATION  OF  THIS  PAGE  (When  Data  Enlar.d) 


REPORT  DOCUMENTATION  PAGE 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


I report  number 


UCB/AM-78-3 ' 


2.  GOVT  ACCESSION  NO 


3.  RECIPIENT’S  CATALOG  NUMBER 


«.  TITLE  (and  Submit) 


n The  Onset  of  Breakup  in  Inviscid 
and  Viscous  Jets 


S.  TYPE  OF  REPORT  6 PERIOD  COVERF.O 

Technical  Report 


«.  PERFORMING  ORG.  REPORT  NUMBER 


7.  AUTMORfJj 

D.  A.  Caulk  and  P.  M.  Naghdi 


• ■ CONTRACT  OR  GRANT  NUMBERfaJ 


N00014-76-C-047V 


».  PERFORMING  ORGANIZATION  NAME  ANO  ADDRESS 

Department  of  Mechanical  Engineering 
University  of  California 
Berkeley,  California  94720 


10.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  0 WORK  UNIT  NUMBERS 


NR  062-536 


II.  CONTROLLING  OFFICE  NAME  ANO  AOORESS 

Fluid  Dynamics  Program 
Office  of  Naval  Research 
Arlington,  Virginia  22217 


IZ.  REPORT  CATE 


July  197S-. 


IS.  NUMBER  OF  PAGES 

26 


U MONITORING  AGENCY  NAME  A AODRESSfU  different  Itom  Controlling  Otllce) 


It.  SECURITY  CLASS,  (ol  CM.  r.por 1) 


ISa.  DECLASSIFICATION/ DOWN  GRADING 
SCHEDULE 


IS.  DISTRIBUTION  STATEMENT  (ol  thlr  Report) 


Approved  for  public  release;  distribution  unlimited 


17.  DISTRIBUTION  STATEMENT  (of  Ida  abstract  ontorod  In  Block  30.  II  dlllotonl  trom  Roport) 


IE.  SUPPLEMENTARY  notes 


IS.  KEY  WORDS  ( Continue  an  torotao  oldo  It  nocoacary  and  Identity  by  block  n umber) 

Instability  of  jets,  incompressible,  inviscid 

and  viscous  fluids,  direct  approach,  surface  tension,  one-dimensional 
theory  of  a directed  fluid  jet,  differential  equations  for  small  motion 
on  uniform  flov,  comparison  with  available  results  from  the  three- 
dimensional  theory. 
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This  paper  is  concerned  with  the  instability  of  inviscid  and  viscous 
jets  utilizing  the  basic  equations  of  the  one-dimensional  direct  theory 
of  a fluid  jet  based  on  the  concept  of  a Cosserat  (cr  a directed)  curve. 
First,  a system  of  differential  equations  is  derived  for  small  motions 
superposed  on  uniform  flow  of  an  inviscid  straight  circular  jet  which 
can  twist  along  its  axis.  Periodic  wave  solutions  are  then  obtained 
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20.  Abstract  (continued) 

for  this  system  of  linear  equations;  and,  with  reference  to  a description 
of  growth  in  the  unstable  mode,  the  comparison  of  the  resulting  dis- 
persion relation  is  found  to  agree  extremely  well  with  the  classical 
(three-dimensional)  results  of  Rayleigh.  Next,  constitutive  equations 
are  obtained  for  a viscous  elliptical  jet  and  these  are  used  to 
discuss  both  the  symmetric  and  the  anti-symmetric  small  disturbances 
in  the  shape  of  the  free  surface  of  a circular-  jet.  Through  a compari- 
son with  available  three-dimensional  numerical  results,  the  solution 
obtained  is  shown  to  be  an  improvement  over  an  existing  approximate 
solution  of  the  problem. 
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